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ZEROS AND FACTORIZATIONS OF QUATERNION 
POLYNOMIALS: THE ALGORITHMIC APPROACH 

VLADIMIR BOLOTNIKOV 


Abstract. It is known that polynomials over quaternions may have spherical 
zeros and isolated left and right zeros. These zeros along with appropriately 
defined multiplicities form the zero structure of a polynomial. In this paper, 
we equivalently describe the zero structure of a polynomial in terms of its left 
and right spherical divisors as well as in terms of left and right indecomposable 
divisors. Several algorithms are proposed to find left/right zeros and left/right 
spherical divisors of a quaternion polynomial, to construct a polynomial with 
prescribed zero structure and more generally, to construct the least left/right 
common multiple of given polynomials. Similar questions are briefly discussed in 
the setting of quaternion formal power series. 


1. Introduction 

Any monic complex polynomial p of positive degree can be represented as 

m 

P=n = lcm (p”l > • • • > Pa k -=z-a k \ a t ec, n fc €N, (1.1) 

k =1 

where a 1 ,... , a m are distinct zeros of p of respective multiplicities n±,..., n m and 
where 1cm means the least common multiple. The set {(a^, Rfc)})ILi may be referred 
to as to the zero structure of p and it is clear that a monic polynomial is uniquely 
recovered from its zero structure by formulas (ED- Since the second representation 
in (11.11) is equivalent to the primary ideal decomposition (p) = fjfcLi(pSfc) f° r the 
ideal (p) generated by p and since this decomposition is unique, the zero structure of p 
can be alternatively defined as the collection of generators of irreducible components 
in the primary ideal decomposition of (p). 

The case of polynomials over a division ring (rather than a field) is not that simple 
eu, nz], m- Here we focus on polynomials over the skew field El of quaternions 

a = x 0 + ixi + jx 2 + kx 3 (x 0 , xi, x 2 , x 3 € M), (1.2) 

where the imaginary units i,j,k commute with M and satisfy i 2 = j 2 = k 2 = 
ijk = —1. For a € El of the form (11.2ft . its real and imaginary parts, the quaternion 
conjugate and the absolute value are defined as Re(a) = xo, Im(a) = ixi+jx 2 + kx 3 , 
a = Re(ct) — Im(a) and |a| 2 = oia = |Re(a)| 2 + |Im(a)| 2 , respectively. 

Two quaternions a and j3 are called equivalent (conjugate to each other) if a = h~ l j3h 
for some nonzero h £ El; in notation, a ~ (5. It turns out (see e.g., 0) that 

a ~ (3 if and only if Re(a) = Re(/3) and |a| = |/3|. 

l 


(1.3) 
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Hence, the conjugacy class of a given a£i form a 2-sphere (of radius |Im(a)| around 
Re (a)) which will be denoted by [a]. A finite ordered collection a = (an,..., a*,) 
will be called a spherical chain (of the length k) if 

«i ~ «2 ~ ~ Oik and OLj+\ for j = 1,..., k — 1. (1.4) 

We let HI [ 2 ] be the ring of polynomials in one formal variable z which commutes with 
quaternionic coefficients. Multiplication in HI[ 2 ] (as well as in H) is not commutative; 
however, any (left or right) ideal in El [ 2 ] is principal. We will use notation 

(p) r ■■= {pq: q £ HI[ 2 ]} and (p) £ := {qp : q £ El [ 2 ]} 

for respectively the right and the left ideal generated by p. The subscript will be 
dropped if the ideal is two-sided; it is not hard to show that any two-sided ideal in 
HI[ 2 ] is generated by a polynomial with real coefficients. 

For any non-constant polynomial / £ HI[ 2 ], the sets Z^{f) and Z r (f) of its left 
and right zeros (see Section 2.1 for definitions) are non-empty and are are contained 
in a finite union of distinct conjugacy classes: 

Z(f) := Zi(f) [J Z T (f) 

Moreover, each conjugacy class Vj either contains exactly one left and one right zero 
of / or Vi C Zi(f) P| Z r (f ); in the latter case, we say that Vj is the spherical zero of 
/. Observe that the set of all polynomials / € HI [ 2 ] having the spherical zero V is 
the two-sided ideal generated by the real polynomial 

X\r(z) = (z — a)(z — a) = z 2 — 2z ■ Re(a) + |aj 2 , a £ V (1.5) 

(the characteristic polynomial of V)\ it follows from characterization (| 1.3 f) that a in 
(11.511 can be replaced by any other element in V. 

Several algorithms for finding left/right roots (by means of complex root-finding 
for a polynomial with real coefficients are known m, m, nu, m- Yet another 
algorithm of this type (Algorithm 13.41 below 1 is our contribution to the topic. 

By the division algorithm, any non-constant monic / £ HI [ 2 ] can be factored as 

/(A) = (z-7i)(z-j 2 )---(z-'yn), 7i> • • • jTat € H, (1.6) 

where 71 £ Zf (f ) and 77 V £ Z r (f). Although it may happen that none of 72 ,..., 7 at-i 
belongs to Z(f ), it is still true that Z(f) C [77]. 

The zero structure of quaternion polynomials is a delicate issue (for related results 
based on various notions of zero multiplicities, we refer to HU, m, [a])- Here we 
propose to characterize the left (right) zero structure of a polynomial in terms of 
its left (right) spherical divisors. To introduce these divisors, we first recall several 
definitions. 

The least right common multiple h = lrcm(/, g) of /, g £ El[ 2 ] is defined as a 
(unique) monic polynomial such that (h) r = (f) r n ( g) r . The least left common 
multiple llcm(/, g) is defined as a (unique) monic polynomial generating the left 
ideal (/)* D (g) t . 

We will say that a (monic) polynomial / is indecomposable if it cannot be repre¬ 
sented as the least right (left) common multiple of its proper left (right) divisors. 
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As we will see in Section 5, a polynomial / is indecomposable if and only if it ad¬ 
mits a unique factorization (II. 6 j) into the product of linear factors and the elements 
( 71 ,... , 7 at) from this factorization form a spherical chain. In other words, there is 
a one-to-one correspondence 

Q! (o\, ■ ■ ■ ■ ®-k) 1 ^ Pol ■ Pa\ * * * Pa^ (f*^) 

between spherical chains and rnonic indecomposable polynomials. Here and in what 
follows we use notation p a (z ) := z — a for a fixed a € HI. 


Theorem 1.1. Given f £ Hfz], let V\,...,V m be distinct conjugacy classes con¬ 
taining zeros of f. Then there exist unique (monic) polynomials F>{ v and D d y 
(i = 1,... ,m) such that 

f = D[ Vi -h i so that Z(D[ Vi ) C V*, Z(h i )f]V i = 9, (1.8) 

f = gi .D f rVi so that Z(D f ryi )cVi, i%)f|VS = 0. (1.9) 

Furthermore, f = Ircm (d £ Vi ,..., F>^ v ) = llcm(ZP) y , ... , D d yZ) or equivalently, 


(f)r — {f)t — 


i= 1 


2=1 


Finally, there exist (unique) integers > 0 and spherical chains on = (ajp,. 
and 5j = (5^ i,..., a^) in V) sncd that 


7-)/ _ y^i 


/«, and I)[ v . = P, 


Ot-i 


■yKi 

*Vi 


for i = 1 ,..., m, 


where P ai and Pg. are indecomposable polynomials defined as in m- 
quently, 


( 1 . 10 ) 

• 5 a i,ki ) 

( 1 . 11 ) 

Conse- 


f = 


Xy™ ' lrcm(P ( 


ct l) 


■ P, 


_ yKi 

— x Vl 


Xyff ' 11 cm (Pj 


CXI ’ ‘ 


.^aj- ( 1 - 12 ) 


We will refer to D'j>y. and Dl v . as to left and right spherical divisors of /. Each 
left (right) spherical divisor contains all information about zeros of / within the cor¬ 
responding conjugacy class and thus, we may define the left (right) zero structure 
of / as the collection of its left (right) spherical divisors. The integer m is called 
the spherical zero multiplicity of the conjugacy class V). The collections {(ac*, aj)}^ 
and {(Ki, 5j)}P 1 of spherical multiplicities and spherical chains from representations 
( 11 . 121 ) also can serve as definitions of left and right zero structures of a given poly¬ 
nomial /. A result related to Theorem 11.11 is the following version of the Primary 
Ideal Decomposition Theorem. 

Theorem 1.2. For any non-constant polynomial f € M[z], there exist left (right) 
relatively prime indecomposable polynomials p \,... ,pm (resp., pi.... ,pu) such that 
f = lcrm(pi,.. . , pm ) = 11cm(pi,... , pm )- Equivalently, 

M M 

(/) r = |>>* and if) l = n<fi>* 

2=1 2=1 


(1.13) 
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In contrast to (11.101) . all ideals in (| 1.13 1) are irreducible, although the polynomials 
Pi and pi in (jl.l3p are not determined from / uniquely. However, the nonuniqueness 
can be described explicitly (see Theorem 15.61 below). 

The objective of this paper is to establish explicit connections between repre¬ 
sentations (11.61) . ( 11 . 101 ) . ( 11 . 121 ) . (11.131) and in particular, to construct explicitly a 
polynomial with the prescribed spherical divisors. One can see from ( 11 . 121 ) . that 
the non-trivial part here is to construct the lrcm and 11 cm of m indecomposable 
polynomials with zeros in rn distinct non-real conjugacy classes. The construction 
has been known for the case where all (but at most one) these polynomials are linear 
(see e.g., ED- The algorithm settling the general case is our next contribution to 
the topic. 

The outline of the paper is as follows. After recalling some known results on 
quaternion polynomials in Section 2, we present the root-finding Algorithm [33] (Sec¬ 
tion 3) and Algorithm 14. 31 which produces the spherical divisor of a given polynomial 
associated with a given conjugacy class. The proof of Theorem 1 1.1 1 is given in Section 
4.1, and two algorithms relating left and right spherical divisors corresponding to the 
same conjugacy class are given in Section 4.2. In Section 5, we give several character¬ 
izations of indecomposable polynomials and prove Theorem 1 1.21 The construction of 
the least common multiple of several indecomposable polynomials is given in Lemma 
15.21 (the case where all polynomials have zeros in the same conjugacy class), Algo¬ 
rithm [674] (where the polynomials have zeros in distinct conjugacy classes). Finally, 
Algorithm 16.61 produces a polynomial with prescribed spherical divisors while the 
least common multiple of given polynomials is the outcome of Algorithm 16.81 In the 
concluding Section 7, we discuss similar questions in the framework of formal power 
series over H. 

2. Background: quaternion polynomials and their zeros 

A straightforward computation verifies that for any a € H and / G H[z], 

f(z) = f ee (a) + (z — a) ■ ( L a f)(z ) = f er {a ) + (. R a f)(z ) • (z - a), (2.1) 

where f ee (a ) and f £r (a ) are respectively, left and right evaluation of / at a given 
by 

m mm 

f e£ {a) = ^2a k f k and f r (a) = ^ f k a k if f(z) = '*T J z?f j , ( 2 . 2 ) 
k =0 k =0 j =0 

and where L a f and R a f are polynomials of degree m — 1 given by 

m— 1 / m—k—1 \ m— 1 / m—k—1 \ 

(W)D=E E aif k+j+ Az\ (R a f)(z)=Y,{ E A+i+iA)A. 

k =0 ^ j =0 ' k=0 ^ j =0 ' 

(2.3) 

Interpreting Ufa:] as a vector space over HI we observe that the mappings f L a f 
and / i-a R a f define respectively the right linear operator L a and the left linear 
operator R a (called in analogy to the complex case, the left and the right backward 
shift, respectively) acting on IHIfz]. 
















5 


2.1. Left and right zeros. An element a. £ HI is called a left (right) zero of / £ BI[z] 
if f ee (a) = 0 (respectively, f er (a ) = 0). We denote by Zt{f ) and Z r (f) the sets of 
all left and all right zeros of /, respectively, and we let Z(f) := Zi(f) U Z r (f). The 
next two equivalences follow from representations ( 12 . 1 : 

/ e£ («) = 0 O / € (p a ) r, r(a) = 0 ^ f£ (p a ) t . (2.4) 

Combining equivalences (12.41) gives the following result (see [4, Lemma 2.6] for de¬ 
tails). 

Lemma 2.1. Let a, (3 £ W be two distinct conjugates: f3 £ [a]\{a}. Then 

(Pa )r ^ (pp )r = {P a )l H {pp)l = (^[a])- (2-5) 

In case the non-real conjugacy class (the 2-sphere) V is a subset of Z(_{f) D Z r (f), 
we will say that V is a spherical zero of /. 


2.2. Polynomial conjugation. The quaternionic conjugation a4aonl can be 
extended to the anti-linear involution / i —> ft on IHIfz] by letting 

m m 

/"(*) = if f( z ) = ^2 zJ fj- ( 2 - 6 ) 

3=0 3=0 

It is not hard to verify that 

//** = / a /, {fgf = {fg)(fg) i = (//“XW*)- (2.7) 

One can see from (|2.21) that if / £ M^], then f ei (a ) = f er {a ) for every a £ H; in 
particular, if / € R[,z], then Z(f) = Zi(f) = Z r (f). Moreover, if f £ M[*], then for 
each a £ H and h ^ 0, we have f(h~ 1 ah) = h~ 1 f(a)h so that Z(f) contains, along 
with each a , the whole conjugacy class [ct]. Since for any f £ H[;z], the polynomial 
//N is real, the set Z(ff^) is the union of finitely many conjugacy classes. The 
following result is essentially due to I. Niven [16]. 

Theorem 2 . 2 . Let deg(/) > 1 and let Z(ff^) = (J Vi be the union of distinct 
conjugacy classes. Then Zi(f) (J Z r (f ) C Z(ff^) and each conjugacy class V) either 
contains exactly one left and one right zero of f or Vi £ Zi(f)f]Z r (f). 

2.3. Zero multiplicities: If we denote by /( fe ) the k- th formal derivative of f £ 
Elfz], then a straightforward verification shows that for any fixed a £ H, 


deg(/) 

/- E 


pi 


(/ (fc) ) e£ («) 


k=0 


k\ 


deg(/) 

E 

k=0 


(. f (k) Y 


a 


k\ 


Pi 


Pa(z) ■= z — a. 


( 2 . 8 ) 


Definition 2.3. Let us say that a £ H is zero of f £ Hfz] of left zero multiplicity 
mi(a] f) = k if / = p^h for some h £ H[z] with h e£ (a) ^ 0 . 

It follows from ()2.8I) that f £ (p„) r if and only if (/0’)) e ^(a) = 0 for j = 0,..., k— 
1 ; therefore, k = mi(a\ f ) can be alternatively defined as the least nonnegative 
integer such that (/( fc )) e *(a) / 0 . 
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The right zero multiplicity m r (a; /) is defined as the integer k such that / = hp ^ 
for some h G Hfz] with h er (a ) 0, or equivalently, as the least nonnegative integer 
such that (/^) er (a) 0 . 

Definition 2.4. Let us say that a conjugacy class V C H is the spherical zero of 
/ £ IHI[xj] of the spherical multiplicity n = m s (V ; /) if 

/0) = Xy(z)g(z) = g(z)Xv(z) (2-9) 

for some g G BI[z] vanishing at at most one point in V. Equivalently, k = m s (V ; /) 
is the integer such that / € (Xy)\(Xy +1 ). 

The local and spherical zero multiplicities are related as follows: for f £ H[z] and 
a conjugacy class V C H, 

m s (V ; /) = min{m^( 7 ;/)} = min{m r ( 7 ;/)}. ( 2 . 10 ) 

7 eV 7 eV 

In fact, as it follows by successive applying Lemma [2.1l to / and its formal derivatives, 
both minimums in (| 2 . 10 |) are attained at all (but at most one) elements in V. 

Remark 2.5. Combining (12.101) with Definition (12.31) of mi and m r we conclude 
that for any two elements a 7 ^ (3 in the conjugacy class V C H, the spherical zero 
multiplicity m s (V; f) equals the least integer k > 0 such that at least one of the ele¬ 
ments (/^) e£ (a) and (f^) ei {fi) (equivalently, one of (f^) er (a) and {f^ K ' > ) er {/3)) 
is non-zero. 


Example 2.6. Let a ~ /3 ~ 7 be in the same conjugacy class V and let 

f(z ) = (z- a) n (z - f3)(z - 7 ) fc . 

(1) # a ^ 7 ^ 7 , then m f (a; /) = n, m r (r,f) = k, m s (V;f) = 0. 

(2) In particular, if (5 7 ^ a = 7 , then mg(a ; /) = n, m r (a; /) = k, m s (V ; /) = 

0 . 

(3) If 7 = P + « + P, then mi(a; f) = n + 1, m r (r,f) = k, m s {V\f) = 1. 

(4) // 7 = a = /3 and k < n + 1, then mi(a; /) = m r (a;/) = n + 1, 
m s (y-,f) = k. 

2.4. Evaluation of products: By (|2.4[) . if g e£ (a) = 0, then ( gf) e£ (a ) = 0 for any 
/ G IHIfz]. On the other hand, since (gf)(z) = z k g(z)fk, we also have 

n 

(gf) e£ (a) = g ei (a) ^^(a)- 1 ™^))* f k = g e * (a)/ e *(<f*(a)" 1 (a)), 

provided g ei (a) 0. Therefore, the left evaluation of the product of two polynomials 
is defined by the formula 

( r ,n e£ try\ - I ffe£ (")' ^ {g ei ( a )~ la a ee ( a )) if g ei (a>)fio, . . 

191 > 0 if g e£ (a) = 0 . 1 j 

Similarly, (gf) er (a) = g er (/ er (a)a/ e 7 , (a) _1 ) -f er (a) if f £r {a) 0 and ( gf) £r {a ) = 

0 if f er (a) = 0 . 
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Remark 2.7. Let 7 = ( 7 i,..., 7 jv) C [ 71 ] = V be a spherical chain. Then the 
polynomial f = p, fl p yt ■ ■ ■ p lN has a unique left zero 71 and a unique right zero 7 ,/v. 

Indeed, since //** = Xy , it follows by Theorem \2.°A that f has no zeros outside 
V. On the other hand, since (a — /3)~ 1 /3(a — f3) = a for any two distinct a, (3 £ V, 
we conclude from ( 12 . 111 ) that for every 7 £ P\{ 7 i},. 

f ei ( 7) = (7 - 7i)(7i - 72)(72 - 73 ) • • • (7iv-i - 7 n) + 0, 
so that Zi(f) = { 71 }. Similarly, one can show that Z r (f ) = { 7 ^}. 

3. A “new” root-finding algorithm for quaternion polynomials 

Several known algorithms for finding roots of a quaternion polynomial by means of 
complex root-finding (see e.g., (16], wi da, m) assume that all complex zeros of 
the real polynomial //*• are known. If x is a real zero of .ff, then x £ Zi(f)r\Z r (f), 
by Theorem 12.21 If a,a € C are complex-conjugate roots of f f' and if f ei (a) = 
f ei (a) = 0, then [a] C Z^(f) D Z r (f), again by Theorem 12.21 The remaining case 
(finding the only left and the only right zero of / in the conjugacy class [a]) is the 
essence of each individual algorithm. Our contribution here is the following. 

Theorem 3.1. Let f £ IHI[^] ? let a and a be complex roots (or any quaternion- 
conjugate roots) of the real polynomial fp, and let us assume that f ei (a) 7 ^ 0 . 
Then the only left root 7 £ and the only right root y r of f in the conjugacy class [a] 
are given by 

7 1 = (a/ e£ («) + a/ e£ (a))(/ e *(a) + / e£ (a)) - \ 

7r = (/ e£ («) - / e *(a)) _1 (a/ e£ (a) - a>f e£ (a)). 

Proof: If a, (3, 7 are any three distinct equivalent quaternions, then 

f e *{ 7 ) =(7 - P)(a - / 3 ) _1 / e£ («) + (a - 7 )(a - py 1 f ei (P), (3.2) 
f er ( 7) =(« - P)~ l f ee -{u )7 - /5(a - P)~ 1 f et (a) 

+ a(a - fd)- 1 f ei (f3) - (a - P)~ x f e£ (P) 7 . (3.3) 

Formula (|3.2I) relating evaluations of the same type was established in [8|; for formula 
m, we refer to [4J Lemma 3.1]. Letting (3 = a simplifies the latter equalities to 

f e£ ( 7 ) =(7 - 7) _1 [(7 - «)/ e *(a) + (7 - a)f ei (a )], (3.4) 

f er ( 7) =(a - a) -1 [/ e£ (a )7 - a/ e£ (a) + af e£ (a ) - f e£ (ay\ ■ (3.5) 

Now we will show that f e£ (a) 7 ^ Zf e£ (a), so that formulas (13.11) make sense. Indeed, 
if f e£ (a) = f e£ (a), then it follows from (13.41) that f e£ ( 7 ) = f e£ (a) for all 7 £ [a]. 
On the other hand, if f e£ (a ) = —f e£ (a), then again by (13.41) . 

f ei { 7) = (7-7) _1 (a-a)/ e "(«)- 

Neither of the two latter cases is possible since f e£ (a ) 7 ^ 0 by the assumption and 
since / does have a left root in [a] by Theorem 12.21 Thus f e£ (a) 7 ^ ±f e£ (a). 
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Let 7 ^ be a left root of /. Then we conclude from (13.411 that 

( 7 1 ~ a)f ei (a) + - a)f e£ (a) = 0 , 

and solving the latter equation for 7 ^ leads us to the first formula in (13.lj) . Similarly, 
for the right root 7 r we conclude from (|3.5jl that 

f ei (a) 7 r - af ei (a) + af e£ (a ) - f e£ (a) 7 r = 0, 

which being solved for 7 r gives the second formula in (13.111 . □ 

Remark 3.2. Formulas (13.11) for unique left and right roots in a given conjugacy 
class V are based on two left evaluations of f at two quaternion-conjugate points in 
V. If f £r (a) / 0, then the formulas for 7 ^ and-y r in terms of right evaluations of f 
are: 

It = (f er (a)a - f er (a)a)(f er (a) - / er (a) _1 , 

7 r = (/ er (a) + / er (a)) _ 1 (/ er (a)a + f £r {a)a). 

The proof is similar to that of Theorem \3.1\ and will be omitted. 

Remark 3.3. If f e£ (a) = 0, then formulas (13.111 take the form 7 ^ = a (as expected) 
and jr = (f e£ (a))~ 1 af e£ (a). Similarly, if f er (a) = 0, then formulas in Remark\3fM 
give 7 r = oi and 7 ^ = / er (a)a(/ er (a)) _1 . 

The root-finding algorithm based on the preceeding analysis follows. Our contri¬ 
bution to the algorithm is part (3d). 

Algorithm 3.4. Given a polynomial f £ H[z] of positive degree, 

(1) Find all complex zeros of the real polynomial ffK 

(2) Each real zero of f f* is a zero of f. 

(3) Evaluate f e£ (a) and f e£ (a ) for each pair {a, a} of complex-conjugate roots 
//“• 

(a) If f e£ (a) = f e£ (a ) = 0, then [a] is the spherical zero of f. 

(b) If f e£ {a) = 0 / f e£ (a), then a £ Z £ (f) and (/^(ajV^/^a) € 

Mf)- 

(c) If f e£ {a) = 0 / f e£ {a), then a £ Z t (f) and ( f e£ (a )) 1 af e£ (a) € 

ZtU)- 

(d) Otherwise, use formulas (13.11) to compute 7 ^ £ [a] 0 Z £ (f ) and 7 r € 
[a] nz r (f). 

Example 3.5. To illustrate Algorithm \3.4\ take f(z ) = z 2 — z (j + 2k) + 2i. Then 

(ff t )( z ') = z 4 + 5z 2 + 4 = (z 2 + 1 )(z 2 + 4) 

and therefore, all roots of f are contained in the union of two spheres V\ = [i] and 
V 2 = [2i]. Since f e£ ( i) = —1 + 2i + 2j — k and f e£ {— i) = — 1 + 2i — 2j + k, the 
formulas (13.111 give the pair of roots in V\: 


= (—2j - 4k)(—2 + 4i) -1 = j, 7r = (-2k + 4j)" 1 (2i + 4) = 0.8k - 0.6j. 



9 


Similarly, evaluating f ei ( i) = —4 + 2i + 4j — 2k and f e£ (— 2i) = —4 + 2i — 4j + 2k, 
we apply formulas (|3.1I) to a = 2i to get another pair of roots 

ll = (—8j - 16k)(—8 + 4i) _1 = 1.6j + 1.2k, 7r = (8j - 4k)" 1 (8 + 16i) = 2k 

in Vi- Note that the obtained roots correspond to two different factorizations of f: 

f(z) = (z - j)(z - 2k) = (z - (1.6j + 1.2k)) (z - (0.8k - 0.6j)) . 


4. Factorizations 


The standard procedure to factorize a (monic) polynomial / £ H[z] of degree 
N > 0 into the product N monic linear factors is the following: starting with 
Qo = f, construct the sequence { 7 j} C H and the sequence of monic polynomials 

{Qj} by 

7j e Zi(Qj- 1), Qj = L lj Q j -\ for j = 1,..., N , (4.1) 

where L a is the left backward shift operator (12.311 . Since deg(Qj) = N — j and since 
Qj- 1 (7j) = 0) it follows that Qj-i = P-yjQj and Qat = 1. We now recursively get 

f = Qo = p lx Q i = p 11 p 12 Q 2 = ... = p 71 p 72 • • • p 7JV - 

The latter factorization is largely non-unique due to non-unique choices of elements 
7 j in (|4.1D . If we will be picking up the elements 7j from a fixed conjugacy class 
V for as long as possible, we will recover the zero structure of / within V. Details 
are furnished in Algorithm 14.11 below. Observe that since //' is a real polynomial, 
m s ([a\]fp) = mi(ar,ffl). 


Algorithm 4.1. Given f € HQ] and a non-real a £ Z(ff^), let k = m s ( [a] ;//**). 
(1) Compute recursively f3 \,..., /?& £ H and Qi,..., Qk £ HQ] by letting Qo = f 
and 

of */ Qf(«) = o, , , 

1 (<*Qf (a) + aQf (a))(Qj e (a) + Qf(a)) Q if Qf (a) ^ 0, 

Qj+i = Lp j Q j for j = 0, 1, (4.3) 

where Lg. is the left backward shift operator defined in (12.31) . T/ien 
(1)/3i,...,/3 fc £ [a]; (2) Z(Q k ) n [a] = 0; (3) / = P^P^ ■ ■ • Pp k Qk- 


(2) Compute recursively f3\, ..., (3k £ H and Qi, ■ ■ ■ ,Qk € HQ] by letting Qo = f 
and 


Pj +1 


ct, 

(Qf(a)-Qf(aO)-f«Qf(a) 


*/ Qf (a) = 0, 

aQf(a)) if Qf(a) Q (1, 


(4.4) 


Qj+i = -Rg. Qj /or j = 0,..., k - 1, 

where /rQ is Pie right backward shift operator defined in (12.311 . Then 
(1) Pi,...,Pk € [a]; (2) Z(Qk) 0 [a] = 0; (3) / = Q fc p^ • • • p^p^. 


(4.5) 
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Proof: According to Theorem l3.ll the element fij+\ defined in (14.2|) is a left zero 
of the polynomial Qj and belongs to [a]. Therefore, for Qj+i defined as in (14.31) . we 
have Qj = pp. +1 Qj + i from which we recursively recover f = Qo = PfaPfo ''' Pp k Qk- 

By properties (|2.7I) . we conclude from the latter representation that //** = XfiQkQl. 

and since k = m s ( [a]; //**), it follows that Qk has no zeros in [a]. This completes 
the proof of the first part of the algorithm. The second part is justified in much the 
same way. □ 

Algorithm 14.11 produces factorizations requested in (11.81) , (11.91) (with Vj = [a], 
D i,[a] = Pp 1 ''' P/ 3 k an d A>] = Pp k ''' PflJ- As was observed in 0, the polyno¬ 
mial Dim (and similarly, Z) r r a i) can be written in a more structured form if some 
consecutive points obtained via (|4.2I) are quaternion-conjugates of each other. If 
(3j + \ = (3j , then py.Py j+l = pp.p-p = X[a]-> an d the latter real polynomial can be 
commuted through all the factors in -Dg [ Q ] to the left. Incorporating this observation 
and Remark 14.21 below, we get a more efficient modification of Algorithm 14.11 

Remark 4.2. For a £ H, a successive application of formulas (12.31) gives 

n—2 n—k—2 / i \ n 

W = W = U +fe , if f{z) = Y J Z k fk- (4.6) 

k=o i=o yi=o J k= o 

Furthermore, if we define the two-terms recursion 

r 0 = 1, ri=2Re(a), r j+ i = rjn - rj-i\a\ 2 for j = 1,2,..., (4.7) 

an inductive argument shows that = ]Cj=o f or a U k > 0. Therefore, 

the formula for L a La in (14.61) depends on Re(a) and |a| rather than a itself and 
therefore, by characterization (11.31) . L a La = LpDp, whenever a ~ f3. It thus makes 
sense to introduce the spherical backward shift operator 

n—2 n—k—2 n 

S [o‘]f = ^2 zk ^2 ri fo+k f( z ) = '52 zk fki ( 4 - 8 ) 

k= 0 i=0 k= 0 

where the real numbers ri are defined in (H3D- 

Algorithm 4.3. Given f € H[z] and a € H, let m s ([a ]; fft) = k. 

(1) Evaluating f^ J ' 1 at a and a, find the least integer k > 0 such that at least one 
of the elements (f^) e£ (a) and (/M) e * (a) is non-zero. 

(2) Compute g = ST, / by K-times application of formula ()4.8I) . If 2 k = k, then 

g has no zeros in [a] and f = = gX\f a y Otherwise, proceed to (3a) and 

(36). 

(3a) Letting Q o = g, use recursive formulas (14.21) . (14.31) k — 2n times to construct 
«!,..., £ H and polynomials Q i,..., Qk- 2 k = P ■ 

Then (or,..., «&_ 2 k) C [a] is a spherical chain, P has no zeros in [a], and 

■I ~ A [ a]Pa 1 Pa 2 ' ' ' Pa k _ 2K P ' 


(4.9) 
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(3b) Letting Q o = g, use recursive formulas (14.41) . (14.51) k — 2 k, times to construct 
5i,..., 5 fc _2 K and polynomials Q i,..., Qk- 2 n = P- 
Then (5i, ... , 5fc_2 K ) C [a] is a spherical chain, P has no zeros in [a] , and 

f = PpS k _ 2K ■■■Pa 2 P Sl *£] • ( 4 - 10 ) 


Proof: The integer k obtained in Step (1) equals m s ([a];/) by Remark 12.51 
Hence the polynomial g = S^f satisfies (12.91) . and m s ([a]-,g ) = 0. On the other 

hand, since //** = X?Agg^ (by (|2.7R ). it follows that m s ([a\; gg^) = k — 2k and we 
may apply Algorithm 14. II (part 1) to the polynomial g to get aq,..., «fc- 2 K €= [a] and 
the polynomial P fc _ 2K with no zeros in [a] such that g = p ai p a2 ■ • • P ak _ 2fi Pk-2 k- If 

a.j -|_i = o.j for some j, then p aj p aj+1 = X\ a \ so that m s ([a];g) > 1 which is a 

contradiction. Therefore aj+\ / oij for all j as desired. Applying the second part 
of Algorithm 14.11 to the polynomial g, we get representation (14.101) . □ 

Remark 4.4. The representation (14.91) is unique in the following sense: if 

/ = X^p n >p a > 2 ■ ■ ■ Pa'G, Z t {G) n [a] = 0, (4.11) 

is another factorization of / with the spherical chain (a' l5 ..., a ' s ) £ [a], then k' = k, 
s = k — 2k, a'j = ctj for j = 1 ,s, and G = P. The representation (14.101) is 
unique in a similar sense. 

Proof: The integers k and k! are both equal to rn s ([a], /) and therefore, k! = k. 
Then we have from (I4.9jl and (14.111) , 

9 Pa.\Pa. 2 ' ’ ’ Pa k — 2K Pk—2K Pa^Pa' 2 ' ' ’ Pa’ s G■ ( 4- I^) 

Therefore gg ** = ~ 4K Pk- 2 kp\ _ 2 k = ^'[a] GG^ and since Pk- 2 n and G have no zeros 

in [a], we conclude that k — 2 k = s. By Remark 12.71 it follows from factorizations 
()4.12[) that oi (and also a \) is a unique left zero of g in V. Therefore, a\ = a\ and 
applying L ai to equalities (14.121) gives 

Pa 2 ’ ' ' Pa s P P° 2 ' ' " Pa' s G- 

Repeating the above argument we subsequently conclude that aj = a) for all j = 
1,..., s and then also P = G. □ 

Corollary 4.5. Any monic polynomial f with all zeros contained in the non-real 
conjugacy class V can be (uniquely) factored either as f = or as 

f = Xy Pai ■p an (otj £ V, a j+1 + aj). (4.13) 

The statement follows from representation (14.91) and Remark 14.31 since the monic 
polynomial P in (14.91) does not have roots and therefore P = 1. 

Example 4.6. To illustrate Algorithm ^. 3[ let us consider the polynomial 

f ( z ) =z^ — (1 + i + j + k)^ + (2 — i + 2j )z^ — (3 + i + 2j + 2k )z 4 

+ (1 — 2i + 4j )z^ — (3 — i + j + k)z 2 + (2j — i )z + i — 1. 
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A straightforward computation shows that 


(//**)(«) = zU ~ 2zU + 8 - U + 27zW - 30z 9 + 50z 8 - 4027 + 55z6 - 30 - 5 

+ 36 z 4 - 12z 3 + 13z 2 -2z + 2 = (z 2 + 1 f(z 2 - 2z + 2). 


VFe see that all zeros of f are contained in the conjugacy classes V\ = [i] and 
V 2 = [1 + i] . V 2 contains isolated zeros of multiplicity one, whereas V\ contains 
zeros of higher multiplicities. Applying Algorithm \4-3\ we first evaluate: f ee ( i) = 


f e£ {~ i) = 0, (/') e *(i) = (/T*(—i) = 0 and 

i) =42i 5 - 30i 4 (l + i + j + k) + 20i 3 (2 - i + 2j) - 12i 2 (3 + i + 2j + 2k) 
+ 6i(l - 2i + 4j) -2(3 -i+j + k) = —8 — 8i — 8j - 24k / 0. 

Thus, f) = 2. The recursion (14.71) takes the form 

ro = 1, ri = 0, r j+ 1 = -rj-i for j = 1,2,... 
and we subsequently get 

S[i]f = z 5 — (1 + i + j + k)z 4 + (1 — i + 2j)z 3 — (2 + j + k )z 2 + (2j — i )z + i — 1, 
9 ■= s [i]f = z 3 ~ (1 + i+j + k)z 2 - (i - 2j)z + i — 1. 

Since m s ([i\, f f^) = 4 and m s ([i\, f) = 2, we proceed to Step 3. We have 


g e£ (i) = 1 + i+j + 3k, g^{- i) = —l + 3i + j —k, (4.14) 


and therefore, by part (3a) of Algorithm ^. 3[ 


<*X =(-i^(i)+i^(-i))(^(i)+^(-i))- 1 = k, 

Q\ =L ai g = L^g = z 2 - (1 + i + ])z + j - k. 


We next compute Qi*( i) = — i + j — 2k, Qx e {— i) = — 2 + i+j and subsequently get 


0L2 =(-iQt £ (i) + \Q\ l (-i))(Qi*(i) + Q\ l (-i)) 1 =j, 
Q2 L-jQx z 1 i. 


The representation (14.91) for f takes the form 


f(z) = (z 2 + l) 2 (z — k)(z — j)P(z), where P(z) = (z — 1 — i). (4.15) 


Applying part (3b) of Algorithmic. 3\ and making use of (14.141) . we get 


ax =(g e *( i) - g ee (—i)) -1 (—ig e *(i) - i^(-i)) = 


2i + j — 2k 
3 



We then compute 



a 2 =(Q?(i) - Qn-ir'HQTO - iQ?(-i)) = 


2i + 26j + 29k 
39 
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and representation (14.1011 for f takes the form 


/(*) = P(z) 



—2i + 26j + 29k \ 
39 ) 



2 i + j — 2k\ 
3 ) 


(z 2 + l ) 2 


(4.16) 


where P(z) = z — 1 — ril + 3 2k • 


4.1. Proof of Theorem 11.11 Let V be a conjugacy class containing zeros of 
a given / E H[z]. If V = {x} where x is a real root of / of multiplicity k, then 
/ = p k h = hp k x are the factorizations requested in (11.81) . (11.91) . In the non-real case, 
representations (11.81) , (11.91) are established by Algorithm 14.31 with 


Df 

U L\ 


= X 


a\Poti ' P 


otk- 2 k 


and D'l , , = X, 


ttrat-2 k 


•P 


ai) 


(4.17) 


proving, therefore, also formulas (jl.lll) . The uniqueness of factorizations ()1.8I) . (|1,9|) 
and (jl.lll) was shown in Remark 14.31 and Corollary 14.51 Since the real polynomial 
f ft has only spherical zeros or isolated real zeros of even multiplicities, it can be 
factored as f ft = riyli where we again let X^ = pf if x £ M. It follows from 
Algorithm 14.31 and formulas (14.171) that deg (D[ v ) = kj so that deg/ = 

Observe that for any right common multiple F of Dp Vi ,..., Dp , the polynomial 


FF^ is a common multiple of relatively prime real polynomials Dpy.Dpy. = Xy, 
(1 < j < m) and therefore, degF > YftjLi kj- Thus, / is a a right common multiple 
of Dg V ,..., Dp Vm (by (11.81) 1 of the minimally possible degree. Therefore, / = 


lrcm(Zl 


/ 


, Dp v ) which proves the first equality in (11.101) . The second equality 


follows similarly. Equalities (11.121) now follow from (11.111) . 


□ 


4.2. The right zero structure versus the left: If / E H[z] is completely factored 
as in (USD, we can construct its spherical divisors using Algorithm 14.31 If / is given 
in the form (14.91) . then its left zero structure is known only within the conjugacy 
class [a]. However, this information is sufficient to recover the right zero structure 
of / within [a]. We recall the backward shift operators L a and Rp defined in (|2.3I) . 

Lemma 4.7. (1) If F E HI[ 2 ] and 7 E HI are such that 2(F) D [ 7 ] = 0, then 

P-yF = Qpp, where /3 = (F ee ftf))~ 1 'yF ei (fy), Q = Rp(p y F). 

Moreover, 2(Q) D [ 7 ] = 0 and 7 = F er ((3)(3(F er (f3)ft l . 

(2) If Q E HI[ 2 ] and (3 E HI are such that 2(Q) D \(3\ = 0, then 

QP /3 = P 7 T, where 7 = Q er (/3)/?(Q er (/3)) _1 , F = LftQpft. 

Moreover, 2(F) D \f3\ = 0 and (3 = (Q et (x))~ 1 "fQ et {x)- 


Proof: The polynomial g = p^F has a unique left zero in [7] (which is 7). By 
Remark 13.31 the unique right zero of g in [7] is given by 

/5 = 5 e£ (7) _ 1 75 e£ (7) = ((7-7)T e£ (7))^ 1 7(7-7)T e ' 6 (7) = (F et -{j))~ l ^F et -ft}). 
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By (12.lj) . g can be factored as g = Qpg, where Q = Rgg. Since (3 is the only right 
zero of g in [7], it follows that Z(Q) n [7] = 0 . Finally, evaluating both parts in 
p 7 F = Qpg at z = /3 on the right gives 

F er (/3)(3 — 'yF £r (f3) = 0 

from which we conclude 7 = F er (/3)/3(F er (/3)) _1 . This completes the proof of the 
first statement of the lemma. The second statement is verified in much the same 
way. □ 

The next algorithm recovers the right spherical divisor of a given polynomial from 
the given left spherical divisor associated with the same conjugacy class. 

Algorithm 4.8. Given f £ BI[z] in the form 

f = D l[a ] P > Where D {[a] = X [a\P ai Pa 2 ' Pa n («j+1 + «?), (4-18) 

let Pq = P and recursively compute 

Sj+i = Pj t -{a n -j)~ l a n -jPf i {a n - j ), P j+1 = R^^p^.Pj) (4.19) 

for j = 0,..., n — 1. Then f can be represented as 

f = PD i[a]’ where D i[a} = P» n ■ ■ ■ P«2Pdn X [a}’ P = P n- ( 4 ‘ 20 ) 


Proof: Based on the first statement in Lemma l4.7l and definitions (14.191) . a simple 
induction argument shows that 

Pa n —j P j = Pj+iPaj and Z(P j+ 1 ) n [a] = 0 for j = 0,..., n - 1. 
Therefore, we get recursively 

P cni ’ ’ ’ Pa n -iPan P — P a\ ' ' ' Pa n -iPctn. P 0 ~ Pai ' ' ' Pa„-i p Pa\ 

~ Pa\ ' ' ' Pa n -2 P 2Po.2Po.\ — ■ ■ ■ — PnPa n Pa n -i ' ' ' Pai' 

Multiplying the last equality by XF (from either side, since £ M[z]) and taking 
into account (14.181) we get (|4.2UI) . Since m s ([a; /]) = k by (14.181) . it also follows that 
5j + i aj for j = 0,..., n — 1. □ 

Example 4.9. Let us consider the polynomial f from Example \f.6\ and its factor¬ 
ization (14.151) . Thus, oi = k, 02 = j and P(z ) = z—l — i. We now apply Algorithm 
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4 - 8 : 


«i = (-j - 1 - i) j(-j - 1 - i) = 


2i + j - 2k 


Pa 2 p o = z 2 - (1 + i+j)«+j -k, Pi = R 5l (p Q2 Po) = Z- 1 - 


i + 2j + 2k 


a 2 = -k — 1 — 


i + 2j + 2k 


-l 


k —k — 1 — 


i + 2j + 2k\ —2i + 26j + 29k 


39 


p ai Pi = z - 


3 + i + 2j + 5k —2 — 2i — j + 3k 


-z + 


P2 = Ra 2 {p a ,Pl) = Z-l~ 


5i + 12k 


'’Ct2\r'a\ - - -^2 

and the factorization f = P2pa 2 Pa 1 ^m coincides with that in (14.161) . as expected. 


We conclude the section with the algorithm recovering D ^ q j from ^ justifi¬ 
cation is based on the second statement in Lemma 14.71 and will be omitted. 

Algorithm 4.10. Given f G El [2] of the form (14.201) . let Po = P and recursively 
compute 

a j+l = P j V ( a n-j) CX n ~jPj Pj+1 = p OLj + \ (PjPa n -j) ( 4 - 21 ) 

for j = 0,..., n — 1. Then f can be represented as in (14.181) with P = P n . 


5. Indecomposable polynomials and irreducible decompositions 

Let us say that a left (right) ideal in H[z] is irreducible if it is not contained 
properly in two distinct ideals of the same type. The generators of irreducible 
ideals, therefore, are the polynomials that cannot be represented as the least right 
(left) common multiple of their proper left (right) divisors. In |17| . such polynomials 
were called indecomposable. In the next theorem, we collect a number of equivalent 
characterizations of indecomposable polynomials. 

Theorem 5.1. Let f G El[ 2 ] be a monic polynomial factored as in (jl. 6 l) : 

f(z) = (z-7i)(z-72)-"(z-7iv), 7i, • • • G El. (5.1) 

The following are equivalent: 

(1) 7 = ( 71 ,..., 7 jv) is a spherical chain. 

(2) 71 is the only left zero of f. 

(3) 7 jv is the only right zero of f. 

(4) (15.11) is a unique factorization of f into the product of linear factors. 

(5) The ideal (f) r is irreducible. 

( 6 ) The ideal (f)i is irreducible. 

Proof: The implication (1) => (2) was verified in Remark l2.71 Let us assume that 
(1) is not in force, i.e., that either 71 9 ^ 7 j or jj+i = 7 j for some j G {1,..., N}. In 
the first case, the conjugacy class [ 7 ^] contains a left zero a of / different from 71 ; 
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in the second case, p 7 p 7 = p 7 p^. = Xy so that / E (TV) and therefore / has 
infinitely many left zeros. This completes the proof of (1) 77 (2). 

Let us assume that (1) holds and let f(z) = (z — j[)(z — 72 ) • • • (z — Y N ) be another 
factorization of / into the product of linear factors. Since both 71 and are left 
zeros of / and since 71 is the only left zero of / by ( 1 ) =7 ( 2 ), it follows that 71 = 7 J. 
We then consider the equality 

f = P'12 '' ‘ Pin ~ Pi'2 ’' ’ P'i'n 

and since still 7 J+ i / 7 j for j = 2,...,1V — 1, we conclude as above that 72 is 
the only left zero of the polynomial L 71 / and that 72 = 72 - We subsequently get 
7 j = 7 '- so that (15.11) is indeed a unique factorization of /. This completes the proof 
of (1)^(4). 

To prove (4) =7 (5), observe that the uniqueness of (15.11) implies that any left 
divisor h of / is of the form h(z ) = (z — 71 )(z — 72 ) • • • (z — j n ) for some re < N (to 
see this, it suffices to compare factorization / = hg with (15.11) 1. Therefore, for any 
two proper left divisors h = p 71 • • • p 7n and h = p 71 • • • p 7fc (re < k < N) of /, we 
have ( h) T D ( h) r = ( h) r ^ (/) r - Therefore, the ideal (f) r is irreducible. 

To prove (5) (2), let us assume that the ideal ( f) r irreducible. If / has zeros 

in more than one conjugacy class, then / is equal to the lrcm of its left spherical 
divisors, by Theorem 11.11 Since in this case, each left spherical divisor of / is a 
proper left divisor of /, it follows that the ideal (/) r is not irreducible which is a 
contradiction. Therefore Z(f) C [a] for some a € H\R. If / = X^ for k > 1, 
then the polynomials g = p 1 ^ and h = p ^ are proper left divisors of / and their 
lrcm equals /. Therefore, (/) r is not irreducible which contradicts the current 
assumption. It now follows from Corollary 14.51 that / is necessarily of the form 
(14.131) . If re = m s ([a ]; /) > 0, then the polynomials 

9 = P ai ■ ■' Pan-iPat 1 and h = Pk 

are proper left divisors of / and their least right common multiple equals / (the 
details are furnished in Lemma [5.41 below: see also Remark 15.51) . We again conclude 
that (f) r is not irreducible which contradicts the current assumption. Therefore, 
re = 0 in representation ()4. 131) and / has a unique left zero. 

We have verified implications (1) 77 (2) =7 (4) =7 (5) =7 (2). Implications (1) 77 
(3) =7 (4) =7 (6) =7 (3) are verified in much the same way. □ 

Notation: In what follows, we will write Vy := {/ E lHI[z] : Z(f) CV} for the 
set of polynomials having all zeros in V, and we will denote by IVy the set of 
(indecomposable) polynomials having one left and one right zero in V. 

Lemma 5.2. Let g, h E IVy (deg(g) = re > k = deg (h)) be given in the form 


9 Pa 1 ' ' ' Pa n > 


h ~ Phi ’ ’ ’ PPk 


(5.2) 
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where a = (aq,..., a n ) and (3 = (j3\,... , (3k) are two spherical chains from the 
conjugacy class V. If g and h are left coprime (i.e., if a± ^ (3\), then 

Xy, if n = k, 

XyPonPai •••Pa n _ fc > if n > k. 

Similarly, if g and h are right coprime (i.e., if a n ^ (3k), then 


f := lrcm(< 7 , h) = 


(5.3) 


/ := 11cm (g,h) = 


XyP 


Xy, if n = k, 

■■■Pa,,, if n> k. 


(5.4) 


ak+i^ak+2 


Proof: Since / is a right common multiple of g and h we have from (15.211 

f = Pai--Pa n -iPa n P = Ppi m "Pl3 k - 1 Pl3 k Q for some p,q£U[z]. (5.5) 

We then use (12.lip to evaluate the latter representations at a\ and at (3\ on the left: 

f e£ (/3i) = {(3\ - on)(dii - a 2 ) ■ ■ ■ (a„_i - a n )p ee (a n ) = 0 , 
f et -(ai) = 0 = (ai - ft) (ft - (3 2 ) ■ ■ ■ (ft_i - Pk)<f l (Pk)- 

Since aq ft, ctj+i 7 ^ aj and (3j + 1 / / 3j , the latter equalities imply p ei (a n ) = 
q ei ((3 k ) = 0. By (|2.4D . p = Pa n Pi and q = P^ k Qi fo r some p±,qi G H[z]. Substituting 
the latter factorizations into (15.51) gives 

f Pai ' ' ' Pa n —iPa n Poi n Pl Pa i ' " " Pa n _iPld-V 

= p 81 ■ ■ ■ Ph-iP 0 kPhk qi = Phi'" Ppk-i^v- 
Applying the spherical backward shift Sy (14.81) to the latter equalities gives 

ft f Pa 1 ’ ' ' Pfil ' ' ' Pftk— i^l" 

Repeating the preceding argument k — 1 more times we get polynomials pi,... ,p k 
and qi ,..., q k such that pj = Pa n _/Pj+\ and qj = Pa k Qj+i, and eventual equalities 

Syf =Pk = Qk (if n = k) or Syf = p ai ■ ■ ■ Pa n _ k Pk = qk (if n > k). (5.6) 

For / to be a right common multiple of the minimally possible degree, it is necessary 
and sufficient that Pk = 1, and then / is recovered from () 5. 6 [i as in (15.31) . The formula 
(15.41) is justified quite similarly. □ 

Corollary 5.3. Let g\,, g m € IVy he such that degfti) > deg(g 2 ) > ... > 
deg (g m ). If gi, ■ ■ ■ ,g m are pairwise left (right) coprime, then lrcmfti, g 2 ,... ,g m ) = 
lrcm {gi,g 2 ) (respectively, llcm(g l ,g 2 ,... ,g m ) = llcm(g 1 ,g 2 )). 

Proof: Let deg (gj) = dj. By Lemma 15.21 / := lrcrftgq, < 7 2 ) = Xyh for some 

h € XVy with deg (h) = ft — d 2 . Since g^g^ = Xy , it follows that Xy is a right 
common multiple of g 2 ,..., g m . Thus, / is a right common multiple of g 2 , ■ ■ ■, g m and 
the least right common multiple of g\ and g 2 . Therefore, / = lrcm(gi,g 2 ,... ,g m ), 
which proves the first statement. The second statement is verified similarly. □ 

By Corollary 14.51 any monic polynomial / G Vy can be uniquely represented as 
the product / = Xyp of a polynomial p G XVy and a power of the characteristic 
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polynomial Xy. Formula (15.31) (along with Theorem 15. II) tells us that the lrcm(g, h ) 
of two relatively prime polynomials g,h G ZVy (deg (g) > deg(h)) is of exactly the 
same form, where k = deg (h) and p is the (unique) left divisor of g of degree 

deg {p) = deg(g) - deg (h) = deg(/) - 2deg(h). 

We thus arrive at the following result. 

Lemma 5.4. Any polynomial f G Vy with m s (V;f) = k> 1 can be represented as 
f = lrcm(#,h), g, he ZVy, Z £ (g) Z £ (h), (5.7) 

and for any such representation, deg (h) = k and deg(g) = deg(/) — k. Moreover, 

(1) If f = Xy, then (15.71) holds for any g,h G ZVy with deg (h) = deg (g) = k. 

(2) If f = Xyp with p = p ai ... p an G ZVy, then all pairs ( g , h) giving rise to 
representation m are characterized by the properties 

(a) deg(h) = k and Zi(h) ^ {«i}; 

(b) g = pq for some q G ZVy such that deg (q) = k and Zi(g) / {a n }. 

Remark 5.5. In Lemma \5-4\ one can choose h = and g = p k a for any fixed 
a € V in case (1) and h = p^ and q = p k an in case (2). Observe that these 
particular choices were used in the proof of Theorem I5.il 

We now formulate a more detailed version of Theorem P Recall that noncom- 
mutative polynomials p\,... ,p n are said to be left (right) relatively prime if each 
polynomial p^ has no common left (right) zeros with the least right (left) common 

f 

multiple of all other polynomials. We also recall that the spherical divisors D) y. 
and 75^y. of a polynomial / are unique by Theorem ll.il 

Theorem 5.6. Let f G lHI[z] be a monic polynomial with spherical zeros V \,..., V m 
and isolated zeros contained in conjugacy classes V m+ 1 ,... ,V n . There exist two sets 
II = {pi}i(j( n and II := {pi}'fiX n of relatively prime indecomposable polynomials such 
that 

m-fin m+n 

if) r = Pi {Pi)r and (f) £ = P (Pi} £ . (5.8) 

i= 1 Z —1 

Representations (15.81) are unique in the following sense: 

(1) For each j = m + 1,... ,n, the set II (resp., li) contains exactly one polyno¬ 
mial from ZVy j (for each j = m + 1,... ,n) which is equal to D^y. (resp., 

D tv P 

(2) For each j = 1,... ,m, the set II (resp., li) contains exactly two polynomials 
from ZVyj, and the least right (left) common multiple of these polynomials 
is equal to Dj, y (resp., D^ v .). 
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Proof: The existence of representations ([5.8ft follows from Theorem 11.11 and 
Lemma 15.41 Let us assume that 


M 

(/)r = P] ( Pi)r 
2—1 


(5.9) 


for a relatively prime collection II = {pj : 1 < i < M} of indecomposable poly¬ 
nomials. In particular, each polynomial pt £ II is relatively prime with the Ircm 
of all other polynomials in II. By Corollary 15.31 it follows that II contains at most 
two polynomials with zeros in the same conjugacy class. Moreover, if II contains 
two polynomials in IVy^ then Vj is a spherical zero of /, i.e., j £ {1 
Otherwise, Vj contains isolated zeros of /, i.e., j £ [rn + l,...,n}. Therefore, 
M = 2m + (n — m ) = m + n. Let us define Fj £ TV, as the lrcm of two elements 
in II n IV v for j = 1,..., m or as a unique element in II 0 IVv for j > m. Then 
we have from (|5.9D . 

n 

</) r =n<^)r, Fj £ V Vj . 

3 = 1 

As we know from Theorem 11.11 the latter representation implies Fj = for 

j = 1 ,,n. This completes the proof of the part concerning the first representation 
in (15.81) . The dual part is verified in much the same way. □ 


6. Least common multiples 

In the two previous sections, we represented a given polynomial / £ H[z] as the 
least common multiple of its divisors of certain type. Now we address the converse 
problem. 

Problem 6.1. Given a finite collection II = {g\,... ,g m } C EI[z], construct explicitly 
lrcm(gi, ... ,g m ) and llcm^ ,... ,g m ). 

Lemma 15.21 and Corollary 15.31 settled the case where II C IVy consists of pair¬ 
wise coprime polynomials. The next result removes the coprimeness assumption. 
Throughout the section, we will be dealing only with right common multiples and 
consequently, with the left zero structure. The dual statements are analogous and 
will be omitted. 


Lemma 6.2. Given g\,... , g m € IVy such that deg(gi) > deg(g , 2 ) > ■ • • > deg(< 7 m ), 

rv 

n 

let g\ | | Pai ■ Pa\ Pa .2 ''' PoL n Id 
i= 1 

9j=Pjhj, where pj = glcd(^, g x ) for j = 2,...,m. (6.1) 

Then 

lrcm(gi, <72; ■ • • ,9m) = XyP ai Pa 2 ''' Pa n _ k : where k = max^deg(hj). (6.2) 


where 
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deg (pj) 

Proof: Since g\ is indecomposable, its left divisor pj is of the form pj = p a ., 

1=1 

by property (4) in Theorem 15.11 Therefore, for a fixed j E {2,..., m}, we have 

rv 

lrcm(^-, 5 i) =p j • lrcmf hj, JJ p a .J 

i=deg(pj)+l 

rv 

n-deg(hj) 

n 

i=deg(pj)+l 

rv 

n-deg(hj) 

= K°‘ {h,) - II (6.3) 

2=1 

The first equality in the latter calculation follows from (16.11) . the second follows by 
applying Lemma 15.21 to left coprime polynomials hj and /o adeg(p • • • p an and since 

the polynomial Xy is real, the third equality follows from factorization of pj, and 
the last equality follows from Lemma 15.41 (see also Remark 15.51) . We now get (16.21) : 

lrcm( 51 , g 2 ,..., g m ) = lrcm(pi, lrcm (g 2 , gi ),..., lrcm (g m , pi)) 

= lrcm ( 01 , lrcm (p^ s(h ' 2) , 51 ) ,..., lrcm (^p^f hm \ 

= lrcm (si, pi) = 

where the first and the third equalities are self-evident, the second equality holds 
due to (|6.3I) . the fourth equality holds due to the choice (16.21) of k, and the last 
equality holds by Lemma I5l2l applied to left coprime polynomials g\ and . □ 

Observe that if g 2 , ■ ■ ■ , g m are left coprime with g±, then pj = 1 for j = 2 ,m 
and Corollary 15.31 follows from Lemma 16.31 We next consider Problem 16.11 for inde¬ 
composable polynomials having zeros in distinct conjugacy classes. The case where 
all polynomials are linear (II = {p 7l ,..., P 7n }), has been known for a while; see e.g., 
[ 6 j. To handle the general case, we need the following preliminary result. 

Lemma 6.3. Let V be a conjugacy class, let F = p ai ■ ■ ■ p Qfe € XVy and let Q € 
Hfz] be such that Z(Q) n V = 0. Define ai,...,5*, E H and Qq, Qi, ■ ■. ,Qk £ H[z] 
by 

Qq = Q. aj = Q* t _ 1 {aij)- 1 ajQj e _ 1 (a j ), Qj = L Qj ( Qj-iP (6.4) 
for j = 1,..., k. Then Z(Qfi) n V = 0 and 

lrcm(T, Q) = FQ k = Qp 5 l p 5z ■ ■ • p 5fc . (6.5) 
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Proof: Since Z(F ) C V and Z(Q) n V = 0, it follows (see [15, Proposition 4.2] 
for the proof) that deg(lrcm(F, Q)) = deg(i ? ) + deg(Q) = deg(F) + k. Therefore, 
it suffices to find a right common multiple of polynomials F and Q of degree equal 
deg(F) + k. 

Based on the second statement in Lemma 14.71 and definitions (16.41) , an induction 
argument shows that 

PajQj = Qj—iPaj and z (Qj) n V = 0 for j = 1 ,..., k. 

Therefore, we get recursively 


FQk Poi\Pol 2 Pak—iPakQk 
Pol\Pcx 2 ' ' ' ^ak_iQk—lPa^ 

Pa\ ’ ’ ’ Pa^._2 Qk—2Potk—i Pock ' ' ' QoPaiPa2 ’ ' ' Pc*k Q Pa\Pa2 ’ ' ' Pctk > 

from which we conclude that the polynomial FQ f, is a right common multiple of F 
and Q. It is clear that its degree equals deg(F) + k. so that (16.51) follows. □ 

The next algorithm produces the lrcm of n indecomposable polynomials 


Fi = P, 


i rQj t 2 


... p 




e ZV Vl (i = l,...,n) 


with zeros in n distinct conjugacy classes V\,... ,V n Cl. 


( 6 . 6 ) 


Algorithm 6.4. Given polynomials (|6.6I) . 

(1) Let Pi := Fi = p ai l p ai 2 ■ • • p ai ki and let i := 2. 

(2) Let Q o := P% and perform the recursion 

a i+l,j = Qj—i( a i+l,j) a i+l,jQj—i(. a i+l,j)i Qj = L ai+l j (^Qj-lPs i+1 j^j ( d -7) 

for j = 1 ,..., ki. 

(3) Let Pi+i ■= PiQki- If i < n — 1, then let i := i + 1 and go to (2). Ifi = n — 1, 
proceed to (4) 

(4) The polynomial P n is equal to the lrcm(i ? i,..., F n ). 


To justify the algorithm, let us assume that P t = lrcm(i ? i ,.... F t ). Then it 
follows from Lemma 16.31 that 

Pi +1 := PiQki = lrcm(Pj,Ti + i) 

= Ircm(Ircm(Fi,..., F t ), F l+] ) = Ircm(Fi,.. .,F i+1 ). 

Since P\ = F\ = lrcm(Ti), we conclude by induction that Pj = lrcm(i ? i,..., iq) 
for all i = 1 ,..., n. □ 

Example 6.5. We illustrate Algorithm \6.f \ by constructing lrcm(p„,p^) for two 
non-real quaternions a /3. We let P\ = p^, and perform two steps of recursion 
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&• 

Qo = Pi = P<>, Q e fa(fa=P 2 -2(3a + a 2 , 

fa = ((5 2 — 2(3a + a 2 )~ l (3((3 2 — 2/3a + a 2 ), (6.8) 

Qi = L /3 (plpfa) = Z 2 + (P - 2a - fa)z + (3 2 - 2/3a + 2a fa - /3fa + a 2 , 
Q\ l {P) = 3/3 2 - 4f3a + a 2 + 2(a - fa fa, 

fa = (3/3 2 - 4/3q + a 2 + 2(a - P)fafa 1 /3(3/3 2 - 4/3a + a 2 + 2(a - /3)/3i). 
Aow we can write the answer: 

G(z ) := lrcm(p 2 , pj|) = (z - a) 2 (z - #l)(z - £ 2 ). (6.9) 

If we choose in the previous example, a = i and f3 = 1 + j, then formulas (|6.8j) give 

/3i = (-1 - 2i + 2j + 2k)“ 1 (l + j)(-l - 2i + 2j + 2k) = 1-^-, 

fa = (-21 - lOi + 50j + 14k)- 1 (1 + j)(—21 - lOi + 50j + 14k) 

— 1588i + 2645j + 980k 
+ 3237 ’ 

and we get the lrcm of polynomials (z — i) 2 and (z — 1 — j) 2 using formula (16.91) . 


We would like to stress that to ensure mi(/3;G) = 2 for G of the form (16.91) . 
we cannot choose fa = fa; actually, it can be shown that lrcm(p^,/r|) = p 2 p| x 
(with fa as in (16.81) ) if and only if a(5 = j3a in which case we also have fa = fa 
Thus, finding least common multiples of special irreducible polynomials p^f. is the 
computational problem of about the same complexity as the generic one. 

We now present several algorithms based on Algorithm 16.41 The first algo¬ 
rithm constructs a polynomial with prescribed left spherical divisors. The in¬ 
put is the collection of polynomials Fj G Vv :j (j = l,...,m), and the output is 
/ = lrcm(i ? i,..., F n ), which, according to Theorem ll.il is a unique monic polyno¬ 
mial such that 

m 

2(f) c{JVj and D^ v . = Fj for j = 1,... ,m. (6.10) 

j =i 

By Corollary 14.51 each Fj is either of the form Fj = p x 3 . (if Vj = {xj} C M) or, 
otherwise, Fj = Xy.pj, where either pj = 1 or pj G IV Vj ■ 

Algorithm 6 . 6 . Let F\..... F m (Fj G V\j) he arranged so that 

Fj = Xy 3 .pj (1 < j < mi) Fj = Xy, (mi < j < m 2 ); Fj = px] (m 2 <j< m), 

wherepj G IVvj for j = 1,... ,mi. Use Algorithm \6.4\ to compute G = lrcm(pi,... ,p mi 
and write 

( m 2 \ / m 2 

n 

j =1 ' y=ml+l 


( 6 . 11 ) 
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Proof: Algorithm 16.41 applies since pi, ■ ■ ■ ,p mi are indecomposable polynomials 
with zeros in mi distinct conjugacy classes. The rest is clear since px] and Xy. are 

J Vj 

real polynomials and since lrcm(iq,..., F mi ) = (n^i X v‘j) ■ lrcm(pi,... ,p mi ). 

□ 

We next consider the case where all polynomials have zeros in the same (non-real) 
conjugacy class V C H. Without loss of generality we may assume that at most 
one of given polynomials is real (i.e., of the form Xy) - if there are several, then 
removing all of them but the one of highest degree will not affect the lrcm. 

Algorithm 6.7. Given polynomials F\ = Xy and Fj = Xy pj (j = 2 
where pj = P QjA P Qja ■ ■ ■ P ajkj € TPv, 

(1) Construct the polynomials gj\ and gj 2 as follows: 

(a) Let gu = pff and g \2 = pfor any fixed a € V. 

(b) For j = 2,..., m, let g 3 \ = PjPa jik . and g j2 = pff . 

(2) Find f = lrcm^i, gj 2 : 1 < j < m) as in Lemma fOl Then we also have 
f = lrcm(f ? i,..., F m ). 

The statement is obvious since by By Remark 15.51 Fj = lrcm (< 7 ^ 1 , g j2 ) for j = 

1 ,m. Note that in Step 2, we may dismiss all polynomials Fj of degree deg Fj < 

K\ (since F\ is a multiple of each such polynomial). Also, if the polynomials pj € 
TVy are pairwise left coprime, then we may choose a in Step 1(a) so that the 
polynomials g 3 i are all pairwise coprime in which case Step 2 simplifies to 

2'. Pick (any) two polynomials in the set {gji,gj 2 : 1 < j < m} with highest 

degrees and find their least right common multiple using formula (US]). 

We finally present the algorithm that produces the lrcm of any finite collection of 
quaternion polynomials. The construction is based on the fact that the left spherical 
divisor of the lrcm(gi,... ,g m ) associated with a conjugacy class V is equal to the 
lrcm of left spherical divisors of gi,..., g m associated with V. 

Algorithm 6 . 8 . Given polynomials g\..... g rn £ H[z], 

(1) Find all conjugacy classes Vj,..., V n containing at least one left zero of at 
least one polynomial from the set. 

( 2 ) Lise Algorithm \j. 3\ to find all spherical divisors D 9 ^ v of gj for j = 1 ,..., m. 

If Z( gj ) r\V t = 0, let D^y = 1 . 

(3) For each i = 1,... ,n, use Algorithm \6.7\ to construct Fi = lrcm (D^y .,..., Dffy ). 

(4) Use Alaorithm \6.6\ to construct f = lrcm(i ? i,..., F n ). 

Then we also have f = lrcm(gi,..., g m ). 

7. Formal power series over quaternions 

As in the commutative case, certain results concerning quaternion polynomials 
can be extended to formal power series over H (see e.g., m, m, Hi)- We are 
particularly interested in power series for which left and right evaluation functionals 
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make sense. We denote by B = {a£i: |a| < 1} the open unit ball in H, and 
introduce the ring 

s OO 

77 = i f(z) = ^2 fj z3 ■ limsu P V\fk\ < 1 
^ j=0 

Observe that for any f(z) = fk% k i n Hr and any a £ B, the quaternion series 
o a> ‘ fk and fk ak converge absolutely, so the evaluation formulas (| 2 . 2 I) 

(with m = oo) and therefore, the notions of left and right zeros (within B) make 
sense. Furthermore, the power series L a f and R a f (defined as in (|2.3I) but with 
m = oo) are also in 77. Equalities (12.ID and therefore, equivalences (11.51) hold true 
in 77 as well as evaluation formulas (12.lip . The conjugate power series p is defined 
as in (12.61) (with m = oo) and relations (12.71) hold for f,g £ 77. Moreover, the real 
power series fp £ M[[z]] belongs to 77 and (by the complex uniqueness theorem) 
has countably many spherical zeroes in B of finite multiplicity each. Theorem 12.21 
extends to 77 as follows: if f € 77, then each conjugacy class V C either 

contains exactly one left and one right zero of f or V C Ze(f) H ^r(f)- 

Remark 7.1. Theorem \3.1\ holds for any f £ T-L. 

Indeed, formulas (13.11) rely on representation formulas (|3.2I) and (13.31) . which hold 
true for monomials fjz 3 and therefore, for elements in R, by linearity. Furthermore, 
since m s ([a], f p) is finite for every a £ B, Algorithms 14.11 and 14.31 apply to / £ 77; 
the only modification is that Qj and Qj obtained via recursions (I4.2I) - (I4.5I) are power 
series from 77 rather than polynomials. Algorithm 14. II recovers in a constructive way 
the following result from [7j. 

Proposition 7.2. Given f £ 77, let V he a conjugacy class such that m s (V, fp) = k. 
There exist (unique) integer k > 0, spherical chains a = (aq,..., aq._ 2k ) C V and 
a = (ai,..., 5 fc_ 2 /c) C V, and power series P,P £ 77 having no zeros in V such 
that 

./ — ^-V Pen pQ 2 ' ' ' Pdk-2K^ — PPdk-2*. ' ' ’ PoL 2 pai^V • 


Although Algorithms 14.81 and 14.101 relating left an right zero structures apply to 
formal power series, they might be efficient only if evaluations f ei and f er admit 
closed formulas. 

Remark 7.3. Proposition |7.4| allows us to introduce spherical divisors of a given 
f £ 77 (which are polynomials uniquely determined from f). Constructing an f £ 77 
subject to conditions () 6 . 10 D (i.e., with finitely many prescribed spherical divisors 
Dg V ) reduces to Algorithm \ 6 . 61 since any such f is of the form f = GH where 

G = lrcm(Z7g y.,..., D’j) ) and H £ 77 has no zeros in B. The latter construction 
does not seem particularly interesting unless we desire to construct an f with some 
extra properties; then the relevant question is how to choose H to achieve this. In 
the concluding section we consider a problem of this sort. 
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7.1 Finite Blaschke products with prescribed zero structure. We denote by 


H 2 = \ h(z) = ^2z J hj : ||/i||| 2 := ^ | hj\ 2 < oo 


3=0 


3=0 


the space of elements in BIffz]] with square summable coefficients and observe that 
H 2 cH. An element / e H 2 is called bi-inner if the equalities 


h||H2 = ||h|| H2 = ||h-/||H2 for all he H 2 


(7.1) 


hold. If / satisfies only the first (only the second) equality in (17.11) . it is called 
left-inner ( right-inner ). It is easily seen that / is left-inner if and only if /' is 
right-inner. 

For each fixed a£l, the power series 


ka(.z) = ^2<x k z k belongs to H 2 and ||k 0 


|2 _ _ 1 


(7.2) 


k =0 


“IlH 2 “ 1—|q| 2 ' 

Letting T a>7 := 1 — (a + a )7 + |a| 2 7 2 and observing equalities 

T a ,7 ' ( ^ = 1 - 7 « and ( ^ a k 'y k V T Qi7 = 1 - 07 (| 7 | < 1 ), 

' k =0 ' ' k =0 2 

we come up with evaluation formulas 

K l (7) = T " 7 • (1 - 70 ), k® r ( 7 ) = (1 - 07 ) • T “ 7 for all 7 GB, (7.3) 

which imply in particular, that k a has no zeros in B. Since k Q ,(z) • (1 — za) = 1, the 
power series k a is the formal inverse of the polynomial 1 — za. Hence, the power 
series 

OO 

b Q (z) := p a (z) • ka(^) = (z — a) ■ ^ a k z k = —a + (1 — \a\ 2 )z • k^(^) (7.4) 

k =0 

can be viewed as the quaternionic analog of the Blaschke factor • It is seen from 
the second representation in ()7.4j) that Zff b Q ) = Z r (b a ) = {a}. In analogy to the 
complex case, b a gives rise to two automorphisms 7 i->- b® £ ( 7 ) and 7 i->- b® r ( 7 ) of 
the closed unit ball B (see m, 0 ). The latter can be derived from the evaluation 
formulas 

b a ( 7 ) = ^a ,7 ' (7(a 2 + 1) - (7 2 + l) a )j Kf(l) = ((« 2 + !)7 - a(7 2 + 1)) • 

which, in turn, follow from (17.31) (with a instead of a) and (17.11) . 

Proposition 7.4. For each a e B, the power series b Q is bi-inner. 

Proof: The left-inner property b Q has been observed in [2j. Our justification is 
similar to that in [2] but does not rely on the Hilbert space structure of H 2 . We first 
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verify that ||b a • /i || H 2 = II^IIh 2 for every h of the form h(z) = d + cz k . Since 

k -1 

b a (z) ■ h(z) = — ad + (1 — |a| 2 ) • ^ oP^dzP + ((1 — |a| 2 )a k ~ 1 d — ac)z k 

3 =1 

OO 

+ (1 — |ct| 2 ) • y] oi^(oi k d + c)^' +fc+1 , 
j=o 

we have, by the definition of the H 2 -norm, 

||b Q • /i||g 2 = |a| 2 |d| 2 + (1 — |a| 2 ) • (1 — |a| 2fc_2 )|d| 2 + |(1 — \a\ : )a k ~ l d — ac \ 2 
+ (1 - |a| 2 ) ■ |a fc d + c| 2 

= |d| 2 + |c| 2 + (1 — |a| 2 ) ■ ^2Re(a fc dc) — 2Re(a fc-1 dca)^ . (7.5) 

Applying the well-known quaternionic equality Re(a6) = Re(6a) to a = a and 
b = a k ~ 1 dc, we conclude from (17.51) that ||b Q • /i||^ 2 = |c| 2 + |d| 2 = ||/i||^ 2 for 
h(z) = d + cz k . Since multiplication by z preserves the H 2 -norm, it follows that 
||b Q • h || H 2 = \\h \\ R 2 holds for all h of the form h(z) = dz t + cz k and therefore, 
for all h £ EI[z]. By the limit argument, the equality holds for every h £ H 2 and 
therefore b a is left-inner. Then b^ = b„ is also left-inner, so that b a is right-inner 
and therefore, bi-inner. □ 

It is clear that the power series (following [2], we will call it a finite Blaschke 
product ) 

B{z) = b ai (z)h a2 (z) ■ ■ ■ b Qm (z)4>, (ai,..., a m £ B, \<f>\ = 1) (7.6) 

is bi-inner. Its zeros are contained in the union of conjugacy classes [ai],..., [a m ] and 
we can use Algorithm 14.11 to construct the spherical divisors of B corresponding to 
each class. Now we will consider the question related to Algorithm 16.61 to construct 
a bi-inner power series with the prescribed (say, left) zero structure. By Remark [7.41 
in case only finitely many conjugacy classes are involved, the latter question reduces 
to the following: given a polynomial 

Ct Pai Pa 2 • • • Pa m > ^1 1 ■ ■ ■ 5 ^ (^*^0 

find a power series H € H 2 with no zeros in B so that f = G ■ H be inner. 

In case atcxj = ajat for i, j = 1 ,..., m (i.e., or, ..., a m belong to the same two- 
dimensional subspace of H), we can take H := k 3l k 32 ... k Sm . Some reduction is 
possible if G has a real zero x (i.e., G = p x ■ G) or a spherical zero [a] (G = 7f[ a ] • G). 
In both cases, it suffices to find H so that G ■ H is inner and then let H = k^ • H 
(in the first case) or H = k a • k„ • H (in the second case). It thus turns out that the 
general problem reduces to the one where G does not have real or spherical zeros. 
The case where a \,..., a m are pairwise non-equivalent (i.e., has m simple roots) has 
been handled in [2], The next theorem settles the general case. 

Theorem 7.5. Given G as in 0Z2D- there exist fit,... ,/3 m (fii £ [ai]) so that the 
power series G ■ kg x • k^ 2 ... k^ m is a finite Blaschke product. 
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Proof: We will assume (without loss of generality) that ai ,... ,a m are all non- 
real and prove the statement by induction. For rri = 1. we let (3\ = ah and the 
statement follows by Proposition 17.41 Assume that we have found 7 + <4 £ [a^] 
(i = 1 ,..., m — 1 ) so that 

Pa 1 ’ Pa 2 ’ ' ' Pam -1 ‘ ' kc5 2 ‘ ‘ ‘ k ( 5 m _ 1 b 71 ■ b 72 • • • b^ m _ 1 tp. (7-8) 

Since evaluation functional at real points is multiplicative, evaluating both sides in 
(|7.8p at zero gives 

= (7172 • ■ ■7m-l) _1 (aia2 • • • 

and since 7 i £ [a*], it follows that \ip\ = 1. Given equality (17.811 and given a m , we 
want to find £ [a*] (i = 1,... ,m) and 7 m £ [a m ] such that 

Pax ■ Pa 2 ■ ■ ■ Pa m -1 ‘ Pam ' ‘ k/3 2 ‘ ‘ ' k/3 m ~ 4l ' ^72 ■ ■ ' b-y m _j • b 7m ?/), (7.9) 

where the unimodular factor ^ is given by ^ = 7 “ 1 </?a m . Multiplying both parts in 
(17.811 by </? _1 b 7rr 4 on the right and comparing the resulting equality with (17. 91) . we 
see that (17.9H is equivalent to 

Pa\ ' ' ' Pa m -1 ’ k<5i ’ ’ ’ k(5 m _iT 3 — Pai ‘ ’ ’ Pa m -i ' Pa m ’ k /3l ’ ’ ’ k /3m> 

which in turn, is equivalent to 

k 5l •••k 5m _ 1 <^- 1 b 7m V’ = • k/^ • • • \tp m . (7.10) 

The construction of /3 \,... ,/3 m , 7 m and t/’ subject to identity (17.101) is as follows: 
letting a\ = a m and <p\ = 1, we compute a' fc and 4 by the recursive formulas 

o4 +1 = (1 - 4 • 4 • (1 - 4 ^dkcpk)- 1 , (7.11) 

4+1 = 4 ■ (1 - 0fc 1 4</>fca / fc )( 1 - 44( 1 44) _1 ) (7.12) 

for k = 2,..., m — 1. Then we let 

Pk = (1 - <4 4*44) ' 4 X 44 ■ (1 - <4 ^fc^fc) -1 (& = 1, • • ■, rn - 1), (7.13) 

j3 m = a ' m , = 7 m = '444 _1 - (7.14) 

It follows from (17.111) that a' k+1 ~ a' k for all k = 1,..., m — 1 and in particular, 
a' m ~ a' x = a m . By (17.141) . we also have /3 m) 7m £ [a m ]. From (|7.13l) and the 
induction hypothesis we conclude that (4 ~ 4 6 [ctfc] for fc = 1,..., m — 1. By 
(17.121) . | <4 + i| = |4| for A: = 1,..., m — 1 and therefore, \4> m \ = 1. We also conclude 
from (17.131) that 

^mPa'm • k /3m = <44^ = (f^lpb^ = ^ -1 b 7m V- (7-15) 

To show that (I7.10p holds we first verify equalities 


44 — 4+i4+i> 44 = 4>k+iPki 4 + 444 — 4+i(i + 4+i4) (7-16) 
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for k = 1,..., m — 1. Indeed, letting for short A^ = 1 — a' k 4> k l 5k4’k, we have from 

ump-tini, 

4 > k+ lCtfc+l = (pk ' (1 typ ^k^k^k^^k^k = ^kC^k^k^k = ^k^ki 

(fok+lftk — ' (1 fife fok^kC^kj&k ^k^k^k — t ^fc^fcAfcAj, — fok&ki 

4> k+ i(l + a' k+1 /3 k ) = 4>k ■ (1 - 0fc 1 40fcai)(Ak 1 (l + Aka'k^dkfpkAk 1 ) 

= ifftk (1 T k^k^k ^k&k^^k 

= 4>k ~ Sk4>ka' k + (4>k ~ Skfika'^ia'k + a4)<^fc 
= 4>k — 5k4> k ot'k + (a4 + ot' k )5k(t>kAkA k 1 
= ‘t'k ~ ^k^kOl'k + 5k4>k(oi'k + a' k ) = (f)k + Skcftk&'k- 

The consequence of relations (17.161) is the polynomial identity 

(1 - z5k)<t> k (z - a' k ) = </>fc+iO - c/fc+iX 1 - zfik), 

which being multiplied by k,5 fc (z) on the left and by k p k (z) on the right, implied 
<j>kP a ' k k/3 fc = k 5fc 0 fe+ i for k = 1,..., m - 1. 

Combining the latter equalities with (17.151) gives 

P am ■ kft • kfo • ka 3 • • • k p m = (f>\p a ’ i • k (3l • kg 2 • k ; g 3 • • • kg m 

= k SxfoPa'z ■ k/3 2 ' kg 3 • • • k.g m 

= k 5l • k fe 0 3 p a , • kg 3 • • • kg m = ... 

= ky • k ^ 2 • • • k Sm-i^mPa 1 ™ • k g m 

= k«i •••k5 m _i<p _1 b 7m ?/; 

which confirms (17.101) and completes the proof of the theorem. □ 
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